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Abstract 

We study both analytically and numerically the nonlinear stage of the instability of 
one-dimensional solitons in a small vicinity of the transition point from supercritical 
to subcritical bifurcations in the framework of the generalized nonlinear Schrodinger 
equation. It is shown that near the collapsing time the pulse amplitude and its width 
demonstrate the self-similar behavior with a small asymmetry at the pulse tails due to 
self-steepening. This theory is applied to both solitary interfacial deep-water waves and 
envelope water waves with a finite depth and short optical pulses in fibers as well. 
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1. According to the usual definition, solitons are nonlinear localized objects 
propagating uniformly with a constant velocity (see, for example, The soliton 
velocity V represents the main soliton characteristics which defines very often the 
soliton shape, in particular, its amplitude and width. The velocity V takes values 
in the bounded interval with the end points given by the conditions of touching in 
the /c-space of the plane u — (k • V) with the dispersion relation of linear waves 
uj = uj(k); intersection of these surfaces means presence of the Cherenkov resonance 
that prevents the existence of stationary localized entities. At these boundaries the 
soliton velocity reaches its extremal value V cr equal to the maximum (or minimum) 
phase velocity of linear waves and, as a result, solitons there undergo a bifurcation, 
either supercritical or subcritical. While approaching the supercritical bifurcation 
point from below or above the soliton amplitude vanishes smoothly according to 
the same - Landau - law (oc \V — Vcr] 1 / 2 ) as for phase transitions of the second 
kind. The behavior of solitons in this case is completely universal, both for their 
amplitudes and their shapes. As V — > V cr solitons transform into oscillating wave 
trains with the carrying frequency corresponding to the extremal phase velocity of 
linear waves V cr . The shape of the wave train envelope coincides with that for the 
soliton of the standard - cubic - nonlinear Schrodinger equation. The soliton width 
happens to be proportional to \V — V cr | -1 / 2 . 

In the case of subcritical bifurcation, the situation is similar to phase transi- 
tions of the first kind: at the critical velocity the soliton undergoes a jump in 
its amplitude. In this case the corresponding theory can be developed near the 
transition point between subcritical and supercritical bifurcations (in analogy with 
the tri-critical point for phase transitions). In the series of papers 0, O, Q, [5] we 
demonstrated that in this case the soliton behavior can be described by means of 
the generalized nonlinear Schrodinger equation (NLSE) for the envelope ip, which 
in the one-dimensional case reads as follows: 

1^ _ a 2 ^ + ^ xx - V + 4z/3M 2 ^ + 7#IVf + 3C|V| V = 
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where ujq = w(ko) and ko are the carrying frequency and wave number, respec- 
tively, A 2 = (V cr — V)/V cr <C 1, uj'q the second derivative of u{k) taken at k = k . 
Here the four-wave coupling coefficient p is assumed to have additional smallness 
characterizing the proximity to the transition from supercritical to subcritical bi- 
furcations. The transition point is defined from the equation p = 0. For example, 
for interfacial deep-water waves propagating along the interface between two ideal 



fluids in the presence of capillarity 
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k 3 , 
i, - A 2 - A 2 

^ ~ 1 + p A 

where p is the density ratio, A = (1 — p)/(l + p) the Atwood number, A 2 r = 5/16 
and p cr = (21 — 8\/5)/ll, as it was shown in Ref. [2|. For p < p cr , the four- 
wave coupling coefficient p is negative, and the corresponding nonlinearity is of the 
focusing type. In this case, solitary waves near the critical velocity V cr are described 
by the stationary (d/dt = 0) NLSE and undergo a supercritical bifurcation at 
V = V cr [2|. For p > p cr the coupling coefficient changes sign and, as a result, the 
bifurcation becomes subcritical. For water waves in finite depth h the coefficient 
p changes its sign at 8 cr = k^h w 1.363 [6] while uj'q is always negative. Thus 
the nonlinearity belongs to the focusing type for 9{= kh) > 6 cr and respectively 
becomes defocusing in the region 9 < 9 cr p, |7|. In nonlinear optics, as shown in 
3)], a decrease of p ("Kerr" constant) can be provided by the interaction of light 
pulses with acoustic waves (Mandelstamm-Brillouin scattering). 

Because of the smallness of p we keep in Eq. (JTJ) the next order nonlinear terms: 
the gradient term (~ (3) responsible for self-steepening of the pulse (analog of the 
Lifshitz term in phase transitions), the nonlocal term (due to the presence of the 
integral operator k, the Fourier transform of its kernel is equal to \k\) and the six- 
wave nonlinear term with coupling coefficient C. Two additional 4-wave interaction 
terms, both local and nonlocal, appear as the result of the expansion of the four- 
wave matrix element Tk^k^ in powers of the small parameters = ki — k : 

T kl k 2 k 3 k 4 = TT + ^-(^ + K 2 + K 3 + K A ) (1) 
Z71 Z7T 
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The existence of the nonlocal contribution in the expansion is connected with non- 
analytical dependence of the matrix element T in its arguments. For interfacial 
deep-water waves (IW) this non-analyticity originates from the solution of Laplace 
equation for the hydrodynamic potential and its reduction to the moving interface, 
or instance, for water waves (WW) with a finite depth the nonlocal term is absent 
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as well as for electromagnetic waves in nonlinear dielectrics 



3fl because of the 



analyticity of matrix elements with respect to frequencies, which is a consequence 
of causality (see, for example, Refs. 0, 8]). In the latter case the spatial dispersion 
effects are relativistically small and can be neglected. 

For both IW and WW near the transition point, u'qC is positive; moreover 7 is 
also positive for IW, and therefore the corresponding nonlinearities are focusing, 
thus providing the existence of localized solutions. Depending on the sign of a 
there exist two branches of solitons. For IW they were found numerically [4], 5] 
using the Petviashvili scheme . Explicit solutions for both kinds of IW solitons 
can be obtained in the limiting case only when V —>■ V cr . For negative \i these 
are the classical NLS solitons with a sech shape. For the subcritical bifurcation 
at V — Vcr the soliton amplitude remains finite with algebraic decay (~ l/W) & t 
infinity 0,0]. When nonlocal nonlinearity is absent (7 = 0) soliton solutions can be 
found explicitly. For both branches at large A the number of waves N — J \ip\ 2 dx 
approaches from below and above the same value N cr which coincides with the 
number of waves N on the solitons with /1 = 0. This property for solitons in fibers 
means that the energy of optical pulse saturates, tending to the constant value 
with a decrease of the pulse duration. 

On the other hand, all solitons of Eq. (CQ) are stationary points of the Hamil- 
tonian H for fixed number of waves : 5(H + A 2 iV) = 0, where the Hamiltonian in 
dimensionless variables is given by 



H 



dx. (2) 



This allows one to use the Lyapunov theorem in the analysis of their stability. 
Here, for the IW, A ~ (V cr - Vfl 2 \p - p cr \-\ C = 319/1281 and 



/i = sign(p-p cr ), p = 6/^427, 7 = 32/^427, (3) 

and for the WW case (compare with Ref. [?]]) 

H = sign(# cr -9), fin -0.397, C « 0.176. (4) 

As shown in , for N < N cr solitons corresponding to the supercritical 

branch realize the minimum values of the Hamiltonian and therefore they are stable 
in the Lyapunov sense, i.e. stable with respect to not only small perturbations but 
also against finite ones. In particular, the boundedness of H from below can be 
viewed if one considers the scaling transformation ip = [1/ a) 1 ^ 2 il) s {x/ a) retaining 
the number of waves N where ip = ip s {x) is the soliton solution. Under this 
transform H becomes a function of the scaling parameter a: 

H{a) = (--^]^ fm A dx. (5) 
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It is worth noting that the dispersion term and all nonlinear terms in H, except 
/ ^\ip\ A dx, have the same scaling dependence oc a -2 . The latter means that at 
p = Eq. (00) can be related to the critical NLS equation like the two-dimensional 
cubic NLS equation. ^From Eq. (jSJ) it is also seen that for /x < H(a) has a 
minimum corresponding to the soliton. Unlike the supercritical case, the scaling 
transformation for the other soliton branch with fi > gives a maximum of H(a) 
on solitons and unboundedness of H as a —>■ 0. Under the gauge transformation 
i/j = ip s e lx , on the contrary, the Hamiltonian reaches a minimum on soliton 
solutions and consequently the solitons with /j > represent saddle points. This 
indicates a possible instability of solitons for the whole subcritical branch, at least 
with respect to finite perturbations. In this paper we investigate this question 
in more details. The main attention will be paid to the nonlinear stage of the 
instability. This problem, indeed, is not trivial in spite of a closed similarity with 



the critical NLSE. It is worth noting that Eq. ([I]) at u = 7 = C = represents an 



integrable model (the so-called derivative NLSE) lOj and exponentially decaying 
solitons in this model are stable. It is more or less evident also that small 
coefficients 7, C cannot break the stability of solitons. This means that in the 
space of parameters we may expect the existence of a threshold. Above this 
threshold solitons must be unstable and the development of this instability may 
lead to collapse, i.e. the formation of a singularity in finite time. 

2. Consider the Hamiltonian (j2J) written in terms of amplitude r and phase <p 

= re iv ): 



H 



f 2,^4 7 2? 2 ^61 1 ( 1 a 5 



dx, (6) 



where by an appropriate choice of the new dimensionless variables the renormal- 
ized constant C = C + f3 2 can be taken equal to 1. Hence one can see that the 
Hamiltonian takes its minimum value when the last term in ([6]) vanishes, i.e. when 



ip x + Pr 2 = 0. (7) 

Integrating this equation gives an x— dependence for the phase, called chirp in non- 
linear optics. It is interesting to note that the remaining part of the Hamiltonian 
does not contain the phase at all. 

First investigate the local model when 7 = 0. Let the Hamiltonian be negative 



in some region Q : Hq < 0. Then, following Refs. [11|, [12], one can establish that 
due to radiation of small amplitude waves H n < can only decrease, becoming 
more and more negative, but the maximum value of according to the mean 
value theorem, can only increase: 

max |^|4>^M. ( 8 ) 

This process is possible only for Hamiltonians which are unbounded from below. 
In accordance with (jSJ) such a situation is realized when \x > 0. In this case the 
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radiation leads to the appearance of infinitely large amplitudes r. However, it is 
impossible to conclude that the singularity formation develops in finite time. 

For 7 > the estimations on the maximum value of are not as transparent 
as they are for the local case. Instead of (jHJ), it is possible to obtain a similar 
estimate, 

. .a S\H\ 
maxl^l > 



N 

However, it is expressed through the total Hamiltonian H and the total number 
of waves N. Besides, two inequalities must be satisfied: H < and iV < For 
interfacial waves, N 2 ~ 1.39035 > N cr « 1.3521. Thus, the maximum amplitude 
in this case is bounded from below by a conservative quantity and this maximum 
can never disappear during the nonlinear evolution. 

Now we consider the situation where the self-steepening process can be neglected 
(/3 = 0). In this case Eq. flTJ becomes 

It is possible to obtain a criterion of collapse using the virial equation (for details, 
Il3l . Il4j). This equation is written for the positive definite quantity 



sec 
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R = J x 2 \tfj\ 2 dx, (9) 

which, up to the multiplier N, coincides with the mean square size of the distri- 
bution. The second derivative of R with respect to time is defined by the virial 
equation 

R tt = %(H-£J (10) 
Hence, for /x > one can easily obtain the following inequality: 

Ru < SH, 

which yields, after double integration, R < 4Ht 2 + ait + a 2 . Here a 1)2 are 
constants which are obtained from the initial conditions. Hence, it follows that for 



the states with negative Hamiltonian, H < 0, there always exists such moment of 
time to when the positive definite quantity R vanishes. At this moment of time the 
amplitude becomes infinite. Therefore the condition H < represents a sufficient 
criterion of collapse (compare with [ill, Q). However, it is necessary to add that 
this criterion can be improved by the same way as it was done in Refs. 3 for 
the three-dimensional cubic NLS equation. From Eq. ([TO]) one can see that for 
the stationary case (on the soliton solution) H s = | / |f/>,s| 4 <ir, in agreement with 
Eq. ([HI). As we demonstrated before for /i > the soliton realizes a saddle point 
of H for fixed N. It follows from (jSJ) that for small a the Hamiltonian becomes 
unbounded from below, but for a > 1 it decreases (this corresponds to spreading). 
Therefore in order to achieve a blow-up regime the system should pass through 
the energetic barrier equal to H s . Thus, for this case the criterion H < must 
be changed into the sharper criterion: H < H s . This criterion can be obtained 
rigorously using step by step the scheme presented in [ji], Q and therefore we 
skip its derivation. 



3. In order to verify all the theoretical arguments about the formation of collapse 
presented above we performed a numerical integration of the NLSE (JTJ) for fi > 
by using the standard 4th order Runge-Kutta scheme. The initial conditions were 
chosen in the form of solitons but with larger amplitudes than for the stationary 
solitons. The increase in initial amplitude was varied in the interval from 0.1% up 
to 10%. The initial phase was given by means of Eq. (JTj). In all runs with theses 
initial conditions we observed a high increase of the soliton amplitude up to a factor 
14 with a shrinking of its width. In the peak region pulses for both IW and WW 
cases behaved similarly. Near the maximum the pulse peak was almost symmetric: 
anisotropy was not visible. The difference was observed in the asymptotic regions 
far from the pulse core where the pulses had different asymmetries for IW and 
WW because of the opposite sign for f3 (see Eqs. (J3j) , (J4j) ) . For the given values of 
(3 we did not observe the simultaneous formation of two types of singularities with 
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FIG. 1: Initial (solid line) and final (dashed line) at t = 1.18 distributions for |^|, 
interfacial waves, self-similar variables. The soliton amplitude was increased by 1%, 
fj, = 1, A = 1. The ratio between final and initial soliton amplitudes in the physical 
variables is about 11. 



blowing-up amplitudes and sharp gradients as it was demonstrated in the rcent 
numerical experiments for the three-dimensional collapse of short optical pulses 



due to se 
equation 



f-focusing and self-steepening in the framework of the generalized NLS 



17 1 and equations of the Kadomtsev-Petvishvili type 18]. 



In our numerical computations we found that the amplitude and its spatial 
collapsing distribution develop in a self-similar manner. Near the collapse point 
in the equation (with > 0) one can neglect the term proportional to fi. In this 
asymptotic regime Eq. ([1]) admits self-similar solutions, 



r(x,t) = (t -t)- 1/4 / 



x 



'.ID 



.(*>- w 

where to is the collapse time. 

To verify that we approached the asymptotic behavior given by Eq. (fill) , we 
normalized at each moment of time the ^—function by the maximum (in x) of its 
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FIG. 2: Initial (solid line) and final (dashed line) at t = 2.7192 distributions for |^|, WW 
solitons, self-similar variables. The soliton amplitude was increased byl%,/i = l,A = l. 
The ratio between final and initial soliton amplitudes in the physical variables is about 
11. 
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FIG. 3: Dependence of 1/ max |^| 4 on time. Interfacial waves. 
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modulus max \if)\ = M and introduced new self similar variables, 

i/>(x,t) = il#(£,r), £ = M 2 (x-x max ), r = In M. (12) 

Here x max is the point corresponding to the maximum of |^|. In comparison with 
those given by Eq. (TTTT) . such new variables are more convenient because they do 
not require the determination of the collapsing time to- 

Fig. 1 and Fig. 2 show typical dependences of | -0 1 as a function of the self-similar 
variable £ at t = (solid line) and at the final time (dashed line) for both the IW 
and WW cases. In both figures one can see a fairly good coincidence between the 
initial soliton distribution and the final one at the the central (collapsing) part of 
the pulse and asymmetry of the pulse at its tails due to self-steepening. The latter 
demonstrates that collapse has a self-similar behavior. The form of the central part 
of the pulse approaches the soliton shape because asymptotically the NLS model 
([1]) tends to the critical NLS system. It is necessary to mention that this has 
been well-known for the classical two-dimensional NLS equation since the paper 



by Fraiman 



id- 



Fig. 3 shows how l/max|?/>| 4 depends on time. This dependence is almost 
linear in the correspondence with the self-similar law ( TTTT) . If the initial amplitudes 
were less than the stationary soliton values, then the distribution would spread in 
time dispersively, in full correspondence with qualitative arguments based on the 
scaling transformations (jSJ). 
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